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Mapping of the magnetic flux density B z (perpendicular to the film plane) for a 
YBa 2 Cu 307 thin-film sample was carried out using a scanning micro-Hall probe. The sheet 
magnetization and sheet current densities were calculated from the B z distributions. From 
the known sheet magnetization, the tangential (B x ,y) and normal components of the flux 
density B were calculated in the vicinity of the film. It was found that the sheet current den- 
sity was mostly determined by 2B x ,y/d, where d is the film thickness. The evolution of flux 
penetration as a function of applied field will be shown. 

PACS numbers: 74.60.Jg, 74.75. +t 


I. Introduction 


Stationary and scanning Hall probes 
have been used for evaluation of high tem- 
perature superconducting (HTS) thin 
films. 1 " 4 The magnetic flux density B z nor- 
mal to the film surface is usually measured. 
Algorithms were developed 4 ' 7 to calculate 
the sheet current density J from the mea- 
sured field above the sample (inverse prob- 
lem). In our previous paper, 4 the inverse 
problem was converted to a magnetostatic 
calculation using 


J = V x M = {j X ,Jy} 

_(dM dM\ (1) 

l dy' dx j 


where the sheet magnetization M is normal 
to the film plane (parallel to z-axis). 

The flux density B z at the point (m,n,z) 
is given by 8 


r 3z 2 — r 2 

M(iJ) J / d bidi 
c. . r J 

Nl,N2 

= i M (i, j)G(m, n, i, j, z) 


B z = m f 

«V= i 4 tc 


( 2 ) 


where r is the distance between the local 


sheet magnetization M(i,j) and the field 
point (m,n,z). N1 and N2 are the total num- 
ber of scanning steps in the x and y direc- 
tions, respectively. The integral in (2) is 
over the area of the grid cell (i,j). Equation 
(2) can be written in matrix notation 4 

B Z =GM. (3) 

where G is a matrix of order of Ni 2 x N 2 2 
and M and B z are column vectors of dimen- 
sions Nl x N2. Equation (3) states that the 
sheet magnetization is uniquely determined 
by the measured B z (x,y). The sheet magne- 
tization M(i,j) allows one to calculate the 
flux density B anywhere around the film in- 
cluding the tangential components B x , y of B, 
which are not obtainable directly from our 
scanning Hall Probe measurements. That 
way one can obtain a full picture of the flux 
penetration. This is a forward problem which 
does not require the time consuming solu- 
tions of inverse matrixes. 

The contributions to J from the gradient 
of B and from the curvature of B were com- 
pared. The purpose of this paper is to illus- 
trate graphically the evolution of flux pene- 
tration into a zero-field-cooled YBa 2 Cu 3 C >7 
(YBCO) thin-film sample. 
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Fig. 1. (a) and (b) B z , (c) and (d) M, and (e) and (f) B x distributions in the remanent state for an ideal film 
with uniform current distributions and for a YBCO thin film (sample Y259), respectively. 





n. Experimental 


The micro-Hall probe was patterned 
from a GaAs quantum-well heterostructure 
thin film (University of Bath). Its active 

area is 25 x 25 (im. The scanning micro- 
Hall probe system is a customized commer- 
cial device developed by Quantum 
Technology Corp. 9 in collaboration with the 
participants from Simon Fraser University. 
Epitaxial quality YBCO thin-film samples 
were prepared on LaAlC >3 substrates by 
pulsed excimer laser ablation from a stoi- 
chiometric target of YBCO. We present the 
data for one of the YBCO films (Y259) in 
this paper. Sample Y259 had a critical tem- 
perature of 90 K, and lateral dimensions = 

1.08 x 1.08 cm, and thickness d ~ 300 nm. 
The Hall probe was maintained at a constant 

height of z ~ 250 fim above the film surface 
while taking the lateral scans. The lateral 
scanning step was 0.3 mm. The experimen- 
tal details can be found in Ref. 4. 

HI. Results and discussion 

Fig. 1 compares the theoretical calcula- 
tions for an ideal film (uniformly distributed 
currents flowing in concentric square paths) 
with the measurement results for the YBCO 
(Y259) thin film. Figs. 1(a) and 1(b) show 
the B z distributions for the ideal film and for 
the YBCO film, respectively. An external 

field H a (jioH a ~ 30 mT) perpendicular to the 
film plane was applied to the YBCO film 
(zero-field-cooled), and then H a was 
switched off. The mapping of B z was car- 
ried out with the film in a remanent state. 
The applied flux B a penetrated fully into the 
film. The ideal film was chosen to have the 
same size as that of sample Y259. The sheet 

current for the ideal sample, j (= J/d) = 2.5 x 

10^ A/cm 2 , was chosen to bring the overall 
dependence of B z close to that of the YBCO 
sample. The distribution of M(x,y) was ob- 
tained by solving the matrix equations (2). 4 
The results are shown in Figs. 1(c) and 1(d). 
Fig. 1(d) demonstrates that the YBCO sam- 
ple was in the saturated state. The contour 
lines of M(x,y) represent the current stream 
lines, and the separation between stream 



X (mm) 


Fig. 2. (a) B z (x) and (b) B x (x) for y = 0 in 
the saturation remanent state; the open cir- 
cles are B x ,z values for the YBCO film and 
the solid lines for the ideal film . 


lines is inversely proportional to the value of 
the sheet current density [Eq. (1)]. 

B x and By at point (m,n,z) are given by 


B 


’x,y 


Nl,N2 

J. % mj) 


■ j 

c. . r J 

i 


(4) 


Figs. 1(e) and 1(f) show that B x (x,y) is dis- 
tributed over two opposing triangles. The 
sample symmetry requires that the compo- 
nent By is distributed in two other triangles. 
In each triangular region, there are either 
B x>z or B y , z components of B . B x and By 
vanish along the film diagonals. 

The distributions of B x , B z , and M for 
the YBCO film are very similar to those of 
the ideal film. However there are noticeable 
deviations of the flux distribution from the 
four-fold symmetry pattern which are caused 
by defects in the YBCO sample, see Figs. 1 
and 2. 
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Fig. 4. Schematic diagram of the integral 
path for the application of Ampere's law. 


Fig. 3. Calculated B x ,z values versus z for 
the ideal film for different grid size of the 
sheet magnetization M(x,y). 

The calculations in Fig. 3 were used to 
demonstrate the dependence of B x and B z on 
the grid size, w, of the sheet magnetization 
M(x,y). B z was calculated at the film center. 
B x was calculated at midway between the 
center and the edge of the film. The calcu- 
lated values of B x are correct if the distance 
above the sample surface z > w. B x de- 
creases for z < w and reaches zero when z — » 
0. This is an artifact of the finite grid size. 
In an ideal sample B x should reach a con- 
stant value by approaching the film surface 
(Ampere’s law). The decrease in B x for z < 
w is due to the change of the sign of B x 
across the mid-plane of the film. The calcu- 
lated values of B z support that view. B z is 
continuous across the film thickness and the 
calculated B z show no noticeable depen- 
dence on w. 

The application of Ampere's Law 

$B-dl = Hoi (5) 

to a path in x-z plane shown in Fig. 4 leads 
to 


AB X AB Z t 

~f + ~fd =fi0J y 


where t and t are the integral paths along the 

x and z directions and AB Z is the difference 
in the B z components which are separated 
laterally by L AB X is the difference between 
the B x components above and below the 


film, see Fig. 4. Since d is comparable to 2X, 

where X is the London penetration depth, j is 
considered uniform over the film thickness. 
(Otherwise, j is the average value over the 
film thickness.) The first term in (6) is due 

to the curvature of B, (BV x B, where B is 

A 

the unit vector along B = B B) and the sec- 
ond term is due to the gradient of B, (VB x 

B ). 10 Since B x and B z are comparable and t 
can be close to d, the flux gradient term is 
much smaller (by approximately three orders 
of magnitude) than the flux curvature term. 
In other words, j is mostly determined by the 
curvature term, 


x,y _ 2Bx,y 

d ~ d 



(7) 


This observation is in agreement with the 
theoretical analysis of Clem. 1 ^ From Eq. (7) 
and from the local sheet current j c = 3 x 10 6 

A/cm 2 corresponding to x ~ -5 mm in Fig. 2, 
it can be shown that the calculated value of 

B x (for z=250 pm) is 13% smaller than that 
obtained from Eq. (7). In fact, the profile of 
the B x (x) curve in Fig. 2(b) is very similar to 
that of jy(x) for the YBCO film. 4 This 
shows clearly that the local current density 
in superconducting thin films determines the 
parallel components of B in the vicinity of 
the film. Fig 2(b) shows that B x does not 
reach the saturated value as in the case of an 
ideal sample. The observed slope of B x in 
Fig. 2(b) is a consequence of the dependence 
of jc on the local B z , see Ref. 4. 
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Fig. 8. B x versus x for the YBCO film (a) in 
different applied fields (b) in a remanent 
state. 

Figs. 5-8 show B z , M, and B x distribu- 
tions for sequentially applied external fields. 
Fig. 9 shows vector field plots of the total 
flux density B in the x-z plane at y = 0 for 
various final states, described in the figure 
caption. The evolution of flux penetration 
and the corresponding current distributions 
can be best viewed by the combination of 
Figs. 5-9. 

Fig. 7(f) shows the non- zero values of By 
over the region where B x should dominate. 
This non-symmetrical flux distribution is 
due to the presence of a defect which is 
centered around x=3, y=0, see Fig. 5(c). 
Fig. 5(c) shows that the sheet current path 
curves around this defect and that causes a 
simultaneous presence of the B x and By 
components in that region. 

Figs. 9(a) and 9(b) show the saturation 
remanent flux distributions for the ideal and 
the YBCO films, respectively. The return 
flux (pointing in -z direction) results in the 
demagnetizing field around the film edge. 
The magnitude of B decreases rapidly away 
from the film surface. 

Figs. 5-9 show that for relatively low 


applied fields (e.g. |ioH a = 2 mT), the in- 
duced supercurrents effectively shield the 
interior of the film from the external flux B a . 
However, the tangential components B x ,y are 
present over the entire film surface. The ex- 
clusion of B z requires screening currents 
across the whole film surface, 4 and the pres- 
ence of Bx,y in vortex-free region is a conse- 
quence of Ampere's law. However note that 
the values of B x ,y in the vortex-free region 
are smaller than those in the vortex-pene- 
trated regions [corresponding to the plateaus 
in B x (x,y) in Fig. (7)] where the supercurrent 
reached its critical value, see Fig. 8. 

Fig. 9(c) shows a noticeable bending of 
the flux lines over the film surface. The 

Lorentz force density F = j x B has two 
components. The tangential component (due 
to B z ) points inwards and is responsible for 
moving the flux lines; this force density is 
balanced by the pinning force density. The 
perpendicular component (due to B x ,y) 
points downwards at the top and upwards at 
the bottom of the film. This part of the 
Lorentz force density causes the tilting of 
vortex lines. 

It is interesting to note that the distortion 
of the applied flux is clearly visible only for 
z < 5 mm, see Fig. 9(c), which is approxi- 
mately equal to the half length of the film 
edge. 

As the applied field progressively in- 
creases, the flux penetrates deeper into the 
film and the region in which flux lines ap- 
preciably deviate from the external flux is 
more and more confined to the center of the 
film, see Figs. 9(d)-9(f); 

Fig. 9(g) shows the flux line patterns in 
the remanent state after the applied field 

|ioH a = 8 mT was switched off. The trapped 
flux shows four hilltops, see Fig. 5(e). It is 
worthwhile to point out that at the hilltops 
(maxima of B z ) both B x , y and J are zero 
[Figs. 8, 9(b), and 9(g)]. The induced cur- 
rents flow in closed loops under each hilltop 
and enclose the hilltops along the sample 
edges [Fig. 6(e)]. This multiply-connected 
current distribution can be reconstructed us- 
ing a linear superposition principle with ap- 
propriate "virgin" states 11 ’ 4 . 
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(a) ideal film: sat. rem. state 
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(b) YBCO film: sat. rem. state 
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Fig. 9. Vector field plots of the flux density B in x-z plane (a) and (b) for the ideal and the YBCO films in the 
saturation remanent state, and for the YBCO film for poH a = (c) 2 mT, (d) 4 mT, (e) 6 mT, (f) 8 mT, and (g) 
0 mT. 
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IV. Conclusion 

The mapping of B z near the supercon- 
ducting film surface is shown here to include 
complete and detailed information of the 
current distributions in the film, and the flux 
density in the film vicinity. The tangential 
components B X) y of the flux density B were 
calculated from the sheet magnetization M 
which was determined from the normal 
component of B, B z . B z was measured by a 
scanning micro-Hall probe. The minimum 
height at which the tangential components 
B x ,y can be evaluated is determined by the 
scan step size. The application of Ampere's 
law indicates that the induced supercurrents 
in the YBCO thin-film sample are nearly 
entirely determined by 2B x ,y/d, which is due 
to the curvature of B. 

The Lorentz force has two components, 
one due to B z , the other due to B x>y . In the 
stationary state, the Lorentz force caused by 
B z is compensated by the restoring force of 
pinning sites and the Lorentz force caused 
by the tangential flux densities B x ,y is re- 
sponsible for the tilting of flux lines. The 
calculation of B in the vicinity of the film 
revealed several interesting features: (a) the 
extent of screening by supercurrents in the 
vertical direction is comparable to that in the 
horizontal direction. With an increasing ap- 
plied field, the flux penetrates deeper into 
the film and the region in which flux lines 
appreciably deviate from the external field is 
more and more confined to the center of the 
films; (b) in the remanent state, the areas 
with a maximum trapped flux density have 
no supercurrents; (c) the tangential compo- 
nents B x and By are mostly present in sepa- 
rate and mutually perpendicular triangular 
regions. A significant presence of B x and By 
in the same region, is a consequence of local 
defects; (d) the tangential flux densities B x , 
By are nearly zero along the film diagonals. 

Acknowledgments 

We thank P.H. Brown and S. Govorkov 
for assistance in testing and assembling the 
measuring system. The assistance of S. 
Bending of University of Bath in providing 
the custom Hall probe is gratefully acknowl- 
edged. 

This work was supported by the Natural 


Sciences and Engineering Research Council 
of Canada, Industry Canada (STP-AIM pro- 
gram), and the BC Ministry of Employment 
and Investment. 

References 

1 M.J. Scharen, A.H. Cardona, J.Z. Sun, L.C. 
Bourne, and J.R. Schrieffer, Jpn. J. Appl. 
Phys., 30, L15 (1991). 

2 P. Briill, D. Kirchgassner, and P. Leiderer, 
Physica C 182, 339(1991). 

3 W. Rauch, H. Behner, and E. Gornik, 
Physica C 201, 179 (1992). 

4 W. Xing, B. Heinrich, Hu Zhou, A. A. Fife, 
and A.R. Cragg, submitted to J. of Appl. 
Phys., 1994. 

5 G.M. Asher, J.T. Williams, C.R. Walters, 
H. Joyce and R.J.A. Paul, IEEE Trans. 
Magn. 18, 540 (1982). 

6 BJ. Roth, N.G. Sepulveda, and J.P. 
Wikswo, J. Appl. Phys. 65, 361 (1989). 

7 P.D. Grant, M.W. Denhoff, W. Xing, P.H. 
Brown, S. Govorkov, I. Entin, J.C. Irwin, B. 
Heinrich, H. Zhou, A.A. Fife, and A.R. 
Cragg, submitted to Physica C, 1994. 

8 J. D. Jackson, Classical Electrodynamics. 
(Wiley, New York, 1975), Chap. 5. 

9 Quantum Technology Corp., 1370 Alpha 
Lake Road, Unit 15, Whistler, B.C. Canada 
VON 1B1, Telephone: (604) 938-0030. 

10 John R. Clem, to be published in the Proc. 
of the 7th int. workshop on critical currents 
in superconductors, Jan. 24-27, 1994, 
Alpbach, Austria. 

11 E.H. Brandt, M. Indenbom, and A. Forkl, 
Europhys. Lett. 20, 735 (1993). 


456 



